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Abstract. We prove that every injective Matrix Product State is the unique ground state of a 
simple hopping theory. We start by studying the low energy spectrum of parent Hamiltonians of 
injective Matrix Product States in a particular long range and system size limit under the validity of 
an asymptotic regime with low particle density. We show that in this limit a natural first quantization 
arises. This allows us to compute a specific type of low energy spectrum. This spectrum depends solely 
on the properties of a quantum channel, i.e. transfer matrix of the ground state, and not on any other 
details of the ground-state. We also review normal quantum channels for which the expression is more 
simplified. 

The construction possibly has some interesting uses for the study of quantum and classical Markov 
processes which we briefly expose. As an application, we revisit the notion of (many-body)-localization 
with our framework. Our calculations revealed that translational invariant Matrix Product States can 
be interpreted as a stationary sea of particles. As a next step rather than starting from some local 
Hamiltonian with random potentials, we consider fluctuations of the local tensors of a continuous one- 
parameter family of Matrix Product States. Localization in 1-dimension, is then understood from a 
simple study of spectral and mixing properties of finite dimensional quantum channels. 
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The general study of strongly interacting quantum many body systems requires a large 
amount of parameters. It appears, however, that a successful description of such systems 
is still possible in a much smaller manifold. Matrix product states (MPS), originating 
from the density matrix re-normalization group (DMRG) algorithm [U [2] , have turned 
out to be one of these particular manifolds used to describe the ground state of these 
systems. Further study of the manifold laid down some ground work for studying lower 
energy properties HE]- Some argumentation can be made for the mathematical and 
physical validity of the methods |S]. However, we argue at the beginning of this work 
that this method is heavily biased by the spectral properties of the transfer matrix of the 
ground state. While every numerical method has its own bias, this opens the question 
as to which condensed matter properties hide inside Quantum Channels, i.e. transfer 
matrix, or trace preserving completely positive operators. 

Recalling Kallen -Lehmann spectral representation in quantum held theory [71 
[8], which connects 2-point correlations with free propagators, it is not completely 
unexpected that indeed the transfer matrix contains some information. 

Apart from exactly solvable and integrable models or Lieb-Robinson bound 
arguments based on the structure of the dispersion relation |B], few models can be 
found with a clear hnite-particle pictures for low excitations. In this work, we propose 
an interesting approach to this problem. We start from a system of size N with a 
dynamic described by a parent Hamiltonian, of an injective matrix product state, with 
hnite range L. We then study the limit of Ai, L ^ oo and possible re normalisation of 
the transfer matrix r(L) under the additional restriction that an asymptotic regime for 
low particle density remains valid. 

Under this limit, a clear particle-picture appears, which is completely related to the 
transfer matrix. Even-more so these excitations do not depend on any other microscopic 
properties of the ground-state. 

In 1992, following the seminal work by Affleck, Kennedy, Lieb and Tasaki [9], 
Fannes, et.al. ng showed that Matrix Product States are the ground states of gapped 
local Hamiltonians. In this work, inspiring from this construction and considering 
an additional limit, we show that not only they are the ground-states of gapped 
Hamiltonians, but also this Hamiltonian describes the hopping bosons in one-dimension. 

This work is divided into two distinct part. In the hrst part ( 2pp ), we elaborate 
on the construction. The details of the main theorem Q can be found the appendix. In 
the second part, we present some applications. First, we briefly describe some possible 
consequences and use for the held of classical and quantum Markovian dynamics (§. 
Finally in the last part ([^, we revisit localization using the framework and philosophy 
of the work presented in the first part. 
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2. Prom Kallen -Lehmann spectral representation to the tangent plane of 
Matrix Product States 


Gunnar Kallen |7] and Harry Lehmann |8] discovered independently that the two point- 
correlation in quantum field theory can be related to the free propagators, 

This is illustrated here for a scalar field theory, where p(/i^) is a positive spectral 
density. For a free scalar field, one can verify that the Fourier transform of the two-point 
correlation function, is inversely proportional to the dispersion relation, 


[ d^xC^^{n\Tcl>{x)(l){y)\n) = - - 

J + le 

It is amazing to see how much information is stored simply inside the correlation function 
of the ground state of the system. 

A first approach for studying low energy excitations in one-dimensional spin lattice 
makes use of the tangent plane of the ground-state MPS in the manifold [H]. For a 
translational invariant MPS |0), 

= ( 1 ) 

i 

A simple proposition for an approximation of the basis of the m-particle subspace are 
the states, 

ni<...<nm 

\(p{B^a\ja}a=l) = J] (2) 

i 

A sub-vector space can be considered, by choosing 5*^*^ = 

A particular approximation of the low energy spectrum is defined from the 
variational problem. 


{i;{B,k}\'4!{B,k}) 


(0|^{H,A:}) = O 


In the case of frustration free Hamiltonians, such as parent Hamiltonians, Hz\4>) = 
0,i/z > 0, both numerator and denominator scale as N. The numerator is then finite 
while the denominator is given by the Fourier transform of the transfer matrix Q, 

- Tr(5t7;[r][5])’ 

Generally maxima of the Fourier Transform of the 2-point correlation function, i.e. 
singular points in the Kallen - Lehmann spectral representation, represent the lowest 
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excited particles. Thus, we can argue that the tangent plane method is biased by the 
spectrum of the Fourier Transform of the Transfer Matrix T. 

However, in the case that this particnlar ansatz is right, spectral properties transfer 
matrix, i.e. Fonrier Transform, yields a certain amonnt of information abont the 
spectrum of the Hamiltonian such as minima of the dispersion relation. 

This particnlar observation leads ns to the starting qnestion of this work. Can we 
extract or make other predictions for the dispersion beyond this one-particle ansatz? Is 
it possible to constrnct a family of models for which the spectrnm can be completely 
extract from the transfer matrix? 

We, indeed, manage to construct a local Hamiltonian for which excitations depend 
pnrely on a qnantnm channel, and of which the m-particle eigenstates can be written 
using the m-particle basis (§. The procedure is based on a limit. It will also be 
interesting to study some consequence when the limit slightly breaks down. 


3. Rewriting Parent Hamiltonians 


Parent Hamiltonians mm of Matrix Product States are dehned as projectors onto 
the kernel of the vector space {Tr(XWj ... Wj+^)}. Uniqueness of the ground state is 
assured if the map X ^ Tr ... A^i+^) is injective. 

In the case of our study, we notice that the following parametrization of the parent 
Hamiltonian, pictnred in hgure 0 , is of nse to ns. Choose C E A4d^, so that 

= E ■ P) 

3 

hi.HL = E (C” 

a,/5 

satisfies the parent Hamiltonian conditions. A simple choice is, 

cg;:g = (4) 










\a)0\ 


We mnst hrst point ont a striking property of parent Hamiltonians which becomes 
apparent from equation ([^ and figure Q. Let us make use of the ansatz equation 
([^ for the choice Ba^ = Aa^B. The action of the parent Hamiltonian, as given 
by eqnations w maps the m-particles wave eqnations onto a subspace spanned by 
1,..., m-particles states. Furthermore, and most importantly, the reader should see that 
the decomposition nnder the action, depends purely on the transfer matrix, F, 


T = ® ( 5 ) 

3 

No other information about the ground state or particnlar Hamiltonian details, is 
encoded inside the dynamical spectrum besides the transfer matrix. We will see that 
we can relate the minima of the dispersion with the spectrum of the transfer matrix in 
the vicinity of the unit circle. 
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Figure 1. Some useful representation for the parent Hamiltonians of a Matrix Product State. The 
choice of the tweaking matrix C is not unique. A possible choice is given by equation For this 
particular choice, the Matrix Product is an eigenvector with eigenvalue 1 of this interaction term. 
Furthermore, they are projectors. However, it is unlikely that some low energy spectrum can be solved 
generically for any of such representation. One notices that all these representations convergence 
towards (g) 1 in the limit of large range and system size. 


Limit, N, L ^ oo Clearly the true excitations generally do not satisfy the ansatz 0. 
However, this is the case in the following limiting procedure. First of all, by a choice of 
gauge we can make sure that the transfer matrix is trace preserving 

with a diagonal fixed point, 

r[l], r[p]=p = diag(A), Tr(p) = l 

In the case of an injective MPS, it can be shown that the fixed point p is unique and 
non-degenerate even more so the spectrum always lies within the unit circle. 

Second, we take the thermodynamic limit L < N ^ oo, while varying the range of each 
local interaction of the parent Hamiltonian. We want a meaningful notion of quasi¬ 
particles. As we add more particles, in the generic case, the energy should increase with 
the number of particles. Only in special cases should bound states appear. A reasonable 
choice is L/N —>■ 0. It should be noted for this particular choice, there is still a notion 
of "local Hamiltonian" under rescaling of the metric with L, as shown in figure Q. 

Finally, we could renormalise the transfer matrix F(L). However, we would like 
the interaction term to forget the correction C, in equation Q. Indeed our choice 
mostly converges as (7 —>■ p~^ 0 1. This is not the case if the second largest eigenvalue 
converges towards the unit circle as 0(1/L). So, we must keep in mind that there are 
some restriction on the re-normalization of the transfer. 
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Figure 2. Since we take the limit L oo, one might argue the physical relevance, as the Hamiltonian 
seems, at first sight, to lose its local property. By blocking sites, represented with red disks, and 
regrouping the interactions (dotted lines) we see however it is not the range of interactions which 
increases, but the dimension of the local Hilbert space (red disks in the second line). While the local 
Hilbert space grows in dimension, the interactions remain local, but unbounded, in the thermodynamic 
limit as L/N —> oo. 


As a side note as long as the spectrum of the renormalized super-operator stays 
within this 1 — 0(l/L)-circle, the martingale method remains valid and the Hamiltonian 
is gapped [TO] . 


4. A Natural First Quantization: Hopping of Virtual Particles 

f.l. One-particle excitations 

We saw earlier that the subspace spanned by 1,..., m-particles states remains stable 
under the dynamic. In the final result of theorem 0 . we prove that they can be 
computed exactly within these subspaces. 

In the introduction, we initially computed the energy using variational methods. As 
we desire to derive the whole spectrum, we seek an alternative approach. The goal is to 
demonstrate that for our limit, the states become exact eigenvectors of the Hamiltonian, 

k) — Aj/v 
(k|k) 

Let us first have a look at the limit without re-normalization of the transfer matrix. 
We can propose the following representation for the one-particle states \'ip{Xa, k}) with 



—)■ 0 ^ we set |k) to be an eigenvector 
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momentum /c, 


\'ip{Xa, k} = Y^ Tr ... ... \i) (6) 

n I 

In Lemma (|^, we demonstrate that such states are indeed eigenstates, with energy 
that can calculated from an effective one-particle Hamiltonian . One of the central 
objects, contained in this effective Hamiltonian, is the Fourier Transform of the Transfer 
matrix 7fc[r], 


Tk\L] — Rp -\- e 


ik 


O Rg O 6 


—ik 


( 7 ) 


2(1 - e*^r) ^ ^ 2(1 - 

This result is very evocative of the Kallen -Lehmann spectral representation, 018]. 
Hence, all information about the spectrum can be derived from the spectral analysis of 

7;[r]. 


The following turns out to be extremely useful. One should notice the invariance 
of the one-particle states with momentum k under the gauge transformation, 

A® ^ - A^'^Ya (8) 

However, this invariance breaks down when additional particles are taken into 
consideration. The results for the many-particle wave-function are examined in the next 
section. We should already announce that in the long range and large system limit, we 
derive the eigenstates using the philosophy of the Bethe ansatz m, Each eigenvector 
of the effective one-particle Hamiltonian, which is related to 7fc[r], is a particle with 
momentum k. These particles can theb be combined in a m-particle wave function. The 
energy is then the sum of the individual energy of each particle, 

H\ij[{Xi, ..., Xk-, fci,..., k^}]) ^ y] E[Xj, kjM{Xi, ..., X,; fci,..., k^}]) (9) 

j 


4-2. Breaking down of the Limit and Bound-States 


One could try to look "beyond" this simple one-particle picture. Generally, there exists 
modes consisting of a superposition of scattering particles. These states consisting 
of more than a single particle, sometimes have a lower energy than the other stable 
individual particles. 

The framework, we have presented here, gives us the perfect opportunity for 
studying this phenomenon and even more so derive how this is reffected in the spectral 
properties of the transfer matrix. 

For this, we need to renormalize the transfer matrix F[L] with the range L. The 
idea would be to let the ground state become "critical". As the spectrum gets closer 
to the boundary of the unit circle, some particles, i.e. eigenvectors of 7fc[F], become 
unstable and splits into other particles. 


Technically, the ansatz breaks down with an error ei, given by equation (A. 14). 


The new excited state consists of the 1-particle state in superposition with its fusion 
options ^ V- 
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For simplicity, and as a sake of illustration, we look at normal unital trace-preserving 
completely positive operators, 

rr* = rT, r*[i] = i, T[Xa] = Xx^Xa 

The eigenvectors are orthonormal and characterized by the structure tensor 

XXi, = E 

C 

Normal unital quantum channels are particularly interesting here as the effective 
Hamiltonian is exactly the Fourier Transform of the Transfer matrix. The eigenvectors 
Xa of the transfer matrix, are finite dimensional representations of the particles with 
wave function \'ip{Xa, k}), 

HmX.,k]) = (2t-2Re^— 

with Xx^ = |Ax,j|- Hence the phase of the eigenvalue of V is interpreted as its 

rest-momentum. 

Continuing our discussion at the beginning of this short section, we let the part 
of the spectrum converge towards the unit circle. The rate of convergence becomes 
clear by comparing {'ip{Xa, kyiHlf^lXa, k}) with the generic one-particle energy, Aa = 
2L - 2 Re .r ; 1 - {^f{Xa: k}\H\f;{Xa, k}), 

a ^a.k ' ■ 

a,13 ’ j=l 

First, we see that Aa > 0, and so decreases the energy. 

When two particles have approximately the same eigenvalues, |Ax„| ~ |Ax;sl, and 
approaches the unit circle as |1 — |Ax„| | ~ 7 (lnT)/I/. The new contribution to the 
energy Aa peaks at k = (fxa ~ 'fxp in the order of Additionally, the contribution 

is proportional to the discrepancy of the one-particle sub-space. Hence for 7 < 1/2, the 
one-particle subspace is unstable under the dynamic. The state Xc, resulting from the 
fusion of Xa and Xj^, has a hnite probability of decaying to these two particles. The 
bound state \'ip{Xc, k}) + Cq\'ip{Xa, k — q; X^, q} emerges. This state has a minimum 
at k = (fxa ~ 4>Xf} ■ 


Example We illustrate this for the transfer matrix with eigenvectors, 1, (T+, o~, 0 ^, with 
respective eigenvalues 1,A+,A_ and A^. We choose A+ = e*'^A = A_, A^ = A^. In the 
long range limit, the low energy spectrum is built by adding particles of the type cr’*', o~ 
or to the wave function. 

Choosing to renormalize the spectrum with the range L, yields a rescaling of the 


errors 
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A(L) = e ei(cr+) = e 2 (cr+) ^ ei(cr^) 27 ln(L)L^ 

7 ln(L) 

For 7 = 1/2, the 1-particle subspace remains stable for \ip{a^,k}), \ip{a~,k}), and we 
get the bound state, 

= |'0{Z,/c}) + J]c(^)|'0{cT+,/c-g;(T“,g} + 4^)|'0{cT“,A:-g;(T+,g} 

4-3. Solving the Low Energy Spectrum in the Long Range Limit 

In 1931, Hans Bethe [12] initiated a new technique which ultimately lead to a resolution 
of the Heisenberg spin-chain. His method, which was later given his name, the coordinate 
Bethe ansatz. In his paper, Bethe introduced n-particles wavefunctions of the form, 

0<Xi<...<Xn^N P^Sn 

where (/){xi, ..., Xn) is a wave-function in first quantization of spins at positions Xi,... ,Xn. 
Such plane-wave decomposition needs to satisfy a consistency condition, 

A{PTj) = S{E'"C ^C''C+i)A{P) 

with Tj the transposition operator of j and j -|- 1. Bethe’s intuition is based on the idea 
that such many-body problems should in the most, yet non-trivial, case only depend on 
the scattering of two particles. Such reduction of many-particle scattering is intertwined 
with the Yang-Baxter equation, 

S{Ui,U 2 )S{u 3 ,Ui)S{u 2 ,U 3 ) = S (u^, U 2 ) S {Ui, U 3 ) S {U 2 , Ui) (10) 

This framework turns out to be the most natural for solving the many-particle 
eigenstates. 

First of all, one should wonder what the dynamic actually does. Lemma’s 0 
offers us an elegant picture. Under the action of an interaction terms, a particle will 
hop to another size at a distance < L. After hopping, this particle is split as a new 
superposition of other particles. This transformation is determined by the Transfer 
matrix and the hopping distance. Secondly, we ask how the interaction affects multiple 
particles simultaneously. When each particles are at the large distance, greater than L, 
from each other, this case never happens. This begs to questions of whether we should 
actually care about such case. We prove in Lemma Q, that, indeed, we should not 
do as such. For such asymptotic regime to valid, it is imperative that the main weight 
of the wave-function consists of particles far apart. This is an additional restriction on 
the allowed density state, meaning that it should be relatively low as it can be found in 
more details in the statement of the Lemma. 
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Which particles should we choose? In the one-particle section, we learned that the 
states are invariant under the gauge transformation (|^. The best choice turns out to 


be, 


= A^^Xa + F, = (Id -e-*'=r*)-i[X„] 

where Ba^ should be interpreted as a particle with momentum ka- 

Generally, it turns out that we cannot find the whole emergent spectrum. However, 
due to the asymptotic regime and the frustration freeness of the Hamiltonian, the 
symmetric states emerge as eigenstates. The coefficients A{P) satisfy, 

= H(P) 

for the transpositions These states represent the particles as hopping bosons. 

As a consequence of the Bethe Ansatz, the energy of a m-particle state is the sum 
of the individual energies, 

E[Xa^,ki; ...;Xa^,km] = ^ 

j 

This implies the scaling of the energy with the particle number, m x 2L. Because in 
physics only differentials of the energy matters, a particle conservation is a consequence 
of this scaling. 

Strangely enough, we can solve the whole energy spectrum for normal unital 
quantum channels. Recalling our past discussion about the one-particle spectrum, 
we saw that there were different particle types for each momentum. The momentum 
was an additional label of the particle type. This is not the case for normal unital 
quantum channels. Thus there are at most D^ — 1 different particles, including their anti¬ 
particles. An interesting consequence is that we are able to solve the low energy spectrum 
exactly. Utilizing the parity symmetry, the spectrum is thus simply any symmetric and 
antisymmetric superposition. Hence, 

A(r„,^p) = ±A(P) 

This complete results are stated formally in Theorem Q and its corolary ([^ with 
their proofs in this appendix, which we summarize modestly as. 

Every Injective Matrix Product State is the vacuum of a Hopping Theory 

5. An Emerging Virtual State Space for Stochastics 

The main subject of this work is centralized around properties of some local Hamiltonian. 
Our tensorial notation of parent Hamiltonians Q showed us that natural local sub¬ 
algebra’s appear in the virtual space of the tensor network. Furthermore, these were 
related only and solely to the transfer matrix. We even showed that the low energy 
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Figure 3. A tensorial representation of the famous 1-dimensional Glauber dynamics often studied for 
Ising. 


spectrum was exactly solvable, thereby interpreting some generators of these snb-matrix- 
algebras as physical particles. This point of view offers particularly interesting research 
directions in the field of stochastics, quantum and classical Markov processes. One conld 
label each of these virtual particles and their position inside the network. As we showed 
earlier, nnder the action of the Hamiltonian, these particles are annihilated at the site of 
the action and new particles are created at the bonndary of the interaction. If it tnrns out 
that such actions happens with a positive probability, then we can talk about a Markov 
process on this particular state space O. We briefly expose some applications for the 
Glauber dynamics, and leave the possibilities open for future work in this direction. 

In 1953, an algorithm was presented for calcnlating a particnlar canonical ensemble. 
This method of computation grew to be known as the Metropolis-Hastings m 
algorithm. This algorithms is based on a Markov chains which takes random samples 
from a probability distribution. A famous use of this algorithms is for studying the 
gibbs state of Ising model, 

Ph = H = ^SjSj+i, Sj e {+1,-1} 

3 

The Matrix Product State representation of the one-dimensional Ising model is given 
by the tensor, 

= Si^aSi,bAa,b, A = ^ ^ ^ 

The most simple Markov chain one can think of is the so-called Glauber dynamics in 
the state space of spin. At each time step, one flips a spin with a probability compnted 
from the configuration of its neighbour. It should be obvious that the generators of 
snch dynamics are local and actually related to a parent-Hamiltonian of some Matrix 
Product State. Therefore, the discussion of this work should offer some interesting new 
perspectives on this very well-known topic. On can check that the local generator. 
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C = Y.j - 1 with given in figure Q, is the 

generator of such Glauber dynamics. In this case, C is given by the following equation, 

C oA^ = 1, {Ao B)ij = AijBij 

with the operation -o being the Hadamard product. In the spirit of the main discussion, 
one can find that Pauli matrix Z is one-and only natural virtual particle. Therefore, 
we get a new state space by using the first quantization language with such particles. 
Writing the state. 


\i), aj e { 0 , 1 } 

i 

similarly to the space of spin configurations a G (+1, —1}^, we have the possibility to 
define a state space which mark the presence of a particle Z at the various sites. As 
promised this new state space, is closed under the dynamics, as one can verify that. 








Tr 


aj — i^aj , 0-^+1 
- 1 ^ 3+1 


..., 1 , ..., 


and. 


^aj — i,aj ,CLj-\-i 


- Tr 
4 ^ 


1 (j ^ 


One sees that the expression A^Z^^A^ is symmetric when aj — 0 and anti-symmetric 
for aj = —1. This implies that there will be no creation of additional particle Z, but 
only a hopping to the left or the right with possible annihilation if another particle is 
already present, 


One dimensional Glauber dynamics are known to be exactly solvable. This small 
calculation yields us the famous relations for the dynamical correlation function, 




— Cn{t^ s) + X{Cn-l{tj s) + Cn+l{tj s)) 


This first quantization, which appears naturally from this notation of the parent 
Hamiltonians 0.0, presents us an interesting new state space. The stochastic 
dynamics, of which the tensor network is a fixed point, acts an error corrector. The 
virtual particles are moved around, created and annihilated a certain rate, while the 
whole converges towards the vacuum, i.e. the Matrix Product State without any 
particles added to the virtual space. 


6. Revisiting Localization 

In 1958, Anderson PD] presented the resnlts of his study of diffusion in random 
lattices. He showed that for a particular form of a Hamiltonian consisting of random 
potentials, when the fluctuations were sufficiently large, diffusion would cease. The 
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Figure 4. |MPS) is a sea consisting of large number particles scattering with one-another. The 
distribution of the momentum has a finite number of peaks corresponding to the phase of the largest 
eigenvalues of the transfer-matrix. Diffusion (red) decreases the variance of the peaks, while localization 
increases it. On the level the eigenvalues of the transfer matrix, this corresponds to respectively pushing 
the eigenvalues towards unit circle and the center. 


result first considered as a strange mathematical trickery, grew to a large field, and 
better understanding of insulation-conductor transition, such as Mott transition 125 

50 years later, a part of the research community focussed their attention on the 
so-called Tensor Networks While much richer phases were discovered beyond 

Landau, the philosophy of local order parameter were shifted towards variation of local 
tensors [23] • Furthermore, Gross-Pitaevski equations have been generalized for such 
states |1|. Time-evolution can be understood as a transport along a Tensor-Network 
manifold. 

Even-though localization is a dynamical property, we show in this work that by 
combining both ideas presented above, the same phenomena is realised on the level of 
the manifold of Matrix Product States. The dynamic is exchanged with variation of the 
local tensors in the manifold, on top of which we consider some random fluctuations. 

A very elegant intuitive interpretation of MPS leads us to this insight. In the 
first part of this work, we showed that the Fourier Transform of the transfer matrix of a 
Matrix Product, contains information about low-energy particles of a local Hamiltonian. 
Therefore, we propose that a translational invariant MPS can be established to be sea 
of scattering particles, from which the momentum-distribution has a finite number of 
peaks, as seen in figure Q, 

IMPS) = |?7r(0)) -F \ {m{k) + m{-k))) -\- \ {m{ki) + m{k 2 ) + m{k 3 ))\ki + k 2 + ks = 0) 

k,m 
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Figure 5. Simple illustration of the different mapping l + l + O^O + l + l^l + O + l. 


This interpretation of MPS becomes more transparent, when trying to explain the 
transition from diffusion to localization as we do in section As shown in [211, by 
studying the z-transform of 2-point-correlation functions we can probe the weight of 
various momenta. In some sense, it means that the phase of an eigenvalue is related 
to the momentum of some set of particles in the sea. We also argued earlier that the 
momentum densities p{k) are related to the eigenvalues of the transfer matrix of the 
MPS, and this was exact for normal transfer matrices. 


7. Set-up 

Consider a continuous one-parameter family of Matrix Product States, {5*[s] = 
ylb)[s](A) + Wb)[s]}, 

\^{B^[s]}) = J^Tr (5*1 [1] ...5*i[A^]) (j+*i ... (12) 

i 

As presented earlier the parameter A G / is the equivalent of the time-evolution, in the 
sense of quasi-adiabatic evolution |2S|- The extra terms in the tensors ITb)[s] are taken 
randomly from some measure (i/u(lTb)[s]). 

We need now to make the claimed equivalence more precise. Let us start with 
the time evolution of the density. The creation of a particle at site 0, followed 
by a measurement the outcome at site n, can be calculated from the correlation 
function. As we see below this can be related to the transfer matrix of the MPS-family 

= (A), 

\i>(0,n)? cc KO„M„) - {Oo)(M„>| 


























Signatures of Lattice Excitations in Quantum Channels 


15 


KWIir"-n|r[M])| 

< ||r”(A) - r~(A)||i 

Anderson studied the time-evolution of the probability distribution, |^(0,n, t)p of a 
particle initially at the origin and evolving under a dynamic described by a Hamiltonian 
H = H + E^W(x). He, [201 EH, showed the amazing result that in one and two 
dimension(s) for finite magnitude of fluctuations, such distribution remains localized in 
the sense that, 

^ \n\ |^(0,n,A)p < C 

n 

As argued above, we can now argue localization when studying the function 5(A, Vb), 


^ \n\ |^(0,n, A)p S(A) 


E 

n 

||r-(A)-r-(A)||i 

E II 

r-(A)-r-(A)||i 


for some measure of magnitude of the fluctuations W. 

The next step, is for us to properly connect A with time t and the diffusion. While 
A maybe defined on a small interval, time always goes to infinity. Therefore, in the same 
spirit as the quasi-adiabatic evolution m, time must be varied much more slowly when 
varying \{t). Additionally, the precise relation \{t) is fixed by the property of diffusion 
of the family of MPS, 

E{X{t),M = 0) < DVi 
The set-up is summarized in the table below. 


Dynamical 

Tensor Network 

W{x) 

1TW[x] 

t 

X{t) : S(A(t),IT = 0) < DVi 

U{t) — exp{itH) 

A»(A) + W«[a;] 

|V’(0,n,t)p 

||r’»(A)-r°°(A)||i 

E„ ||r-(A)-r-(A)||i 


8. Localization 

Finally, we can proceed and describe localization. While proving Anderson’s localization 
is extremely challenging, this mapping onto a Tensor Network-picture gives a fresh 
simplistic view of the problem. The following explanation is illustrated in figures 0. 
and ([^. Let us bring our insight of the translational MPS ( |ll| ) come to play. The reason 
for diffusion in the first place is that the leading terms of the wave function is some wave- 
packet. This is compatible with our picture as this means that the distribution p{k) of 
momentum has sharp peaks around a finite number of momenta. Localization arises, 
when leading terms of the wave-function consists of a large number of different wave- 
packets with different momentum. This is again compatible with the sea-picture m}. A 
translational invariant MPS has zero-momentum. In the case of a sea of localized states. 
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Figure 6. Plot of S(A) vs t for the example Q with N = 100 and W = 1. 

the distribution p[k) becomes more and more peaked around the origin. Therefore no 
particles are scattering as, so they are all localized. 

Let us translate the sea-picture, O-l-l+l, for the transfer-matrix, i.e. l-hO+l. We 
saw in the previous section that the probability distribution |t/^(0, n, t) p is related to the 
spectrum of the transfer matrix T. 

The probability distribution p{k) can in principle be studied through by doing 
some tomography and next studying the 2 ;-transform of the 2-point correlation function. 
Keeping this idea in mind, it is acceptable to see the phase of the eigenvalues of the 
transfer matrix, as the momentum of a set of particles in the MPS-sea. Clearly then 
diffusion arises when the eigenvalues are close to the unit circle. Next, local fluctuations 
increase the variance of the peaks. This is, again, can be understood from the spectrum, 
as this means that the eigenvalues are moved towards the origin. The spectrum of 
generic transfer matrices have often a star-shape. This is mathematically understood for 
completely positive operators, as the eigenvalues on the unit circle are a representation 
of some abelian group This fits perfectly with the interpretation of the MPS 

as a stationary sea. The total momentum is zero, and so the sum of the phase must be 
zero, as it is for any non-trivial one-dimensional representations of abelian groups. In the 
extreme case of localization all eigenvalues are zero except one. remaining eigenvalue, 
which is 1, is the zero-momentum of the sea. 

Let us illustrate this with an example. 


Example 1. Consider the local tensors, 

-|- = ^\/l — A -|- wi/V^ , 

- A + W 1 /V 2 J a 

A^^\\) + = (Va + (j^ + W2I21 


with Wj ~ 


N(0, W) and mutually independent for some finite variance W. Here we take 
A G [1/2,1], and see that AKLT belong to the family for A = 2/3. By choosing the right 
function \{f) : [1, 00 ] —>■ [1/2,1] : t —>■ we can plot S(A, M) in figure (6). A simple 

calculation can verify, 


S(A(t), W = 0)< C"(l + Vi), E{X(t), W=1)<C 
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As promised, we see in figure (|^ that our choice \{t) yields diffusion. Similarly to 
Anderson’s result, a finite magnitude of the fluctuation leads to localization. 

As mentioned earlier, localization is especially interesting for studying mott 
insulators. Hence, we need to be able to extend the result to fermions. This can 
of course be done for Fermionic MPS [H5) by replacing a"''*'' with fermionic 

creation operators. The transfer matrix may vary depending on the local observable, 
r = A° ® ® or r = A° ® A° — A^ ® A^. The rest of the set-up is then the 

same, this t im e however, we need to choose transfer matrix with the largest 5(A,M). 

9. Conclusion 

In this work, we proved that all injective Matrix Product States are the unique 
ground-states of hopping models. This was done by studying the spectrum of parent 
Hamiltonians in the long range and system size limit, and eventual renormalization of 
the transfer operator of the Matrix Product ground state, under the restriction that an 
asymptotic regime remains valid for low particle densities We showed that, in this limit, 
the spectrum depends solely on the properties of a quantum channel, i.e. the transfer 
matrix. The derived formula for the low energy spectrum was very reminiscent of the 
Kallen -Lehmann spectral representation in quantum field theory. Particles types are 
represented by eigenvectors of the fourier transform of the transfer matrix 0 - Higher 
energy wavefunction are constructed by putting these particles in the virtual space of 
the ground-state with different momenta. 

Secondly, we consider a renormalization of the spectrum of the transfer matrix 
with range L. The choice was made in such a way that the ground state becomes 
"critical". It results that the particles derived earlier are not necessarily stable. The 
first quantization breaks down when a particle can be fused from two others, of which 
the respective eigenvalues are close to the unit circle. 

This work introduces an interesting interpretation of translational invariant Matrix 
Product States as nothing but stationary sea of particles scattering with one-another. 
This insight motivated us to understand Anderson localization on the level of Tensor 
Networks. The property is introduced to the Matrix Product State manifold without 
considering any Hamiltonian-dynamic. The local random potentials are exchanged with 
fluctuations of the local tensor of the network. The time-evolution is reinterpreted 
as transport along the manifold which moves the eigenvalues of the transfer towards 
the unit circle. This is shown under proper rescaling of the parameter with t to be 
equivalent to diffnsion. Finally, it is illustrated how finite magnitude of fluctuation kills 
the diffusion similarly to the result proven by Anderson. 
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The technical details of the limit are presented in this section. We start from a 
one-dimensional spin chain of size N with interactions of range L. We derive a few 
convergence rates for L ^ oo for the low particle number eigenstates and their energy 
under a possible choice of re-normalisation of the transfer matrix r(L). 

We need to compute the the action of the Hamiltonian onto a m-particle subspace. 
As, we see a very natural new basis appears in the virtual space. However, this basis 
has to be orthogonalized. We, therefore, will have to get some information about the 
Gram-matrix of each subspace. 

First, we start our discussion with results about the one-particle states. In Lemma Q 
we demonstrate that the particle have a simple finite-dimensional representation closely 
related to a quantum channel. Next, when combining these particles to construct multi¬ 
particle states, it can be assumed that the particles are very far from each other. Since 
the dynamic is frustration free, this implies that the theory is non-interactive, and the 
total energy is the sum of the energies of the individual particles. 

Lemma 2 (One-particle energy). Given a Matrix Product State \(/)), the one-particle 
states with momentum k, \'il!{Xa,k}), given by, 


IHXa, n}) = Tr ... ... \i) (A.l) 

i 

are eigenvectors of the local Hamiltonian H{L), with as choice of for the tweaking 
matrix C the relation 0, 


The values E[Xa, k] are the solutions of the finite-dimensional eigenvalue problem of the 
matrix, 


{h{ab)) - 


Tr (p-^rk[xi]rk[xi 


Tr 


,Xa = J2Ua,bXh, U^U = 1 


The basis vectors Xa the eigenvalues of the Fourier transform of the transfer matrix 


0 . 

TilW] = 


For momentum k 
vacuum. 


0, there is the additional constraint of orthogonality onto the 
Q*p°TkO Qp[Xa] = Qp[.] = [■] — P 


Proof. By our choices given by the equations ([^ and Q, the Hamiltionan H{L) is local 
and frustration free. Let us fix the momentum k. Our first step is to show that the 
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one-particle subspace spanned by states of the type \ip{Xa, k}) is closed under the action 
of H{L), in the limit L ^ oo. Knowing this, we can proceed and find an orthogonal 
basis for this subspace, and compute the effective Hamiltonian. 

A look at the choice of corrections for H[L) to be frustration free, yields that in 
the limit L ^ oo all these corrections should converge, 

C —y p ^ ® 1, L —y oo 


Using the frustration freeness of the Hamiltonia we can bound this first error. 


H\f,{Xa,k}) 


jd<n<j+L 




where the correction term of the interactions terms are given by p~^ ® 1. This 

error typically decays as 0(D^|A|^), with A the second largest eigenvalue of the transfer 
matrix T. 

This simplifications leads to the actions, 

= 2L\^P{Xa,k}) - ‘2\^p{H[Xalk}) 

jJ<n<j^L 

- Y J + L}) Tr (etC{A„, k}[e^]) (A.2) 

with, 

L-1 

k}[.] = Y eikjR,-. O p O Lx^ o 

J = 1 

V‘['1 = ['K‘. Op[J = [-l-pTrll, ix[] = Q,[X.QA\] 


An important super-operator is the Fourier Transform which is here approximated by, 




Akn 


IZp o T '"' e 


^—ikn 


The first line in the equation (A.2), is an effective actions on the finite-dimensional 
subspace spanned by Xa- However, the states |^{Aa,A:}) are not necessarily 
orthonormal. This is sorted out by computing the Gram-matrix of the 1—particle 
subspace, 

G^\Xa,XQ = NA[X:,kMXb,k]) = Tr 


The super-operator Tk\L] is hermitian and has thus orthonormal eigenvectors. Using 
the eigenvectors, we can construct the effective one-particle Hamiltonian presented in 
the statement of the Lemma. One should note that the identity particle Xa = 1 with 
momentum A: = 0 is nothing but the ground-state. However, the particle is not an 
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eigenvector for A: 7 ^ 0 . As all one-particle states need to be orthogonal onto the vacuum 
state, for A: = 0, the additional constraint for A: = 0 has to be imposed, 

(^|tA[Xfc,A:])=Tr(pAfc) = 0 

We should finally discuss the error of the action. 




N 


ATV (aIt^Vkaj) 


, 0, eft, L]\(t,[ea„i, e*, j + L]) 


TV (eL.C{A„, A^Xeftl) TV (ei/{A„, (;}[eft]) 


Typically, this error dies off as 0(D^|A|'^). As all the errors are being reduced 
exponentially with L without any dependence on the system size, the claim of our 
lemma is proven. 

□ 


In the theory of scattering one compares the change of the wavefunction of two 
particles initially infinitely far from each other, with the resulting outcome after infinite 
time when they are again far apart. The total energy of the system, since both particles 
are uncorrelated in both cases, is the sum of the individual energy. This additive feature 
of the energy appears regularly in exactly solvable models or even integrable ones. The 
particles simply fill in some Fermi sea. 

We prove in the next lemma, that this is also exactly what happens for parent 
Hamiltonians in the long range and large system size limit. The intuition is that the 
action of the Hamiltonian onto a m-particle wavefunction, approximately only evaluates 
the energy when the particles are very far apart. Clearly for such asymptotic regime to 
be valid, particles need enough empty space. Thus, there is a restriction on the particle 
number which scales with L. 


Lemma 3 (Asymptotic Regime at Low Densities). For a m-particle state with particles 
satisfying the conditions of Lemma 

\HXe.,K}a=l)= E 


ni<...<nr) 


l^{A„,n„C.,) = 52 Tr (A*!... ... W'*-... \t) 


(A.3) 


the action of the Hamiltonian takes into consideration the individual particles 
independently from the others, 

^ 2mL|tA{5«, + 2\f^{H[Xr], kp,X2, k^, ...; A;^}) 

T •. • T I , ki ,..., X^Yi—x , kyn—i, "AA[A^], A)^}) 

L 


The error scales typically as O I — n + 1) 


]\j--n ]j2(n+l) 
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Proof. The proof is more cumbersome in its notation than complex. As a first 
illustration, let us first look at what happens for 2 particles-basis states. One can 
verify, 

= (4L|^[X„,,fci;X,„A:2]) (A.4) 

- \marM-,nkAXa,]M])) (A.5) 

- A:]) (A.6) 

- E E A;e^, A + L; A„„A]) Tr {eiC{X^,M}[eg\) (A.7) 

- E E A; e„, A, 6/302 + L]) Tr (etA: 2 }[eA) (A.8) 

<A,/5 0<ji<j2,|jl—j2|>^ 

- E e/3,i + L]) Tr (e^, fci; AIM) + e* 

with. 


(A.9) 


C(2){X,,,A:i;A„„A;2}[.] = 


T.[-] = E 6 -*^^rA.] 

J=1 


Abui+fc 2 i 2 ]^ O F'l O Lx o F 2 -h o L; 


E e"''""-'" O O i^X ui-"-Uivx , 

p P^ai P^a2 I 

o<ii<i2<n 


The lines (A.4 A.5) are the sought results for the 2-particles. The line (A.6) resnlts from 


the one-particles being transported due to the interactions and fusing with the other 


particle. Lines (A.7 A.8) are the errors of the one-particle approximation, we argned 


in the previous lemma. The last term (A.9) results from the simultaneous action of 


the interaction on both particles. We have already demonstrated that (A.7- and (A.8) 


decay exponentially with L. Hence, only (A.6) and (A.9) are new. One should notice 


that when normalized, the overlap of mi with m 2 particle states scales as. 


^l+'>^2 - 


|WA'„,«:„}™i|V.{V„X„}2,)| = 0(N— 


This is elaborated in more details in the proof of the theorem. The result is similar for 


the state in (A.9) and taking it to be similar to a one-particle state. Looking specifically 


at each line. There are at most two new 1-particle states in the second line. Each new 
one-particle state is the result of LD^ other terms. The argumentation is similar for 


(A.9). There are D^L{L — 1) contribution to this one-particle state. We omitted the 
error terms e* which come from the local identity operation of the interactions. 


e* = E 2(1^1 - n2\ - m; ^ 2 } 

ni l<n2—ni<L 

Under normalization with {'ip[Xaj^,ki;Xa 2 ,k 2 \\p[Xaj^,ki;Xa 2 ,k 2 \), this error is of order 

0{D^E^/N). 

















Signatures of Lattice Excitations in Quantum Channels 


22 


With this in mind, we can tackle the m-particle case and regroup the different error 
types together, 

“1“... , k\ ,..., kYii—x , LLI^Xyh^ , kyii^'j ) 

2 

— II 5 ^2 5 ^2) ••• ) ^nfi) ^m\ ) H“ ••• “h 5 ^2) ••• ) ^m—1-) ^m—1) || 

+ II |C{^2,l{^l5 ^ 2 }; -^3) ^ 3 ; ••• ; km}) + ••. 

|C{^2 5 ^2) ••• ) ^m—2) km—2) ^2,m—l}k m—l) km };^i,fci;})ll 

+ ••• + II IC{‘^m-l,l{^l) ^2) ••) kjn-l}, Xjnkm}) + ... + \'tp{Xm-2,2{k2, km}',Xi, fci}) || 

+ |||C{‘^m,l{^l) ^2) ..) ^m})|| + 


L 

n 


jj 2 (n+i) ^g];.] 2 is. Under the states 


The operator consists of at most 

\C{X{qi ,..., Qrai)] Xi, ki ]...; Xm 2 i km 2 })i 0116 should understand, 
|C{T;Xi,A:i;...;X„,fcJ) 

^ e%i+...+9^Jni ^ ^ n^,+i]) Tr (e^ ^[e;^]) 

ni+L<...<nm^+i,a,g 


with bounded super-operator X. Again the error e* is shown to be, 

n\ ni<...<nm,\ni-nm.\<L 


This yields an error of this order 0{mD‘^^/N'^) Similarly to our illustration the 
overlap of the states xp and (p are bounded as some power of N. Summing over all 
contributions and taking into account the normalization of the m-particle state yields 
the claim. □ 


This Lemma mathematically implies that the Hamiltonian decouples for this basis 
set. In Quantum Mechanics, the stationary physical states, i.e. eigenvectors of some 
Hamiltonian, are orthonormal. Thus, our final step is to orthonormalize this basis. This 
seems quite challenging at first. We could of course use a Gram-Schmidt procedure to 
orthogonalize each m-particle subspace. However, this does not solve the orthogonality 
between each subspace. The representation used previous lemma seems too weak to 
finalize the result. One should for each m-particle eigenvector take a linear superposition 
of not only m-particle states but also l,2,...,m-l. Solving this problem seems very hard 
at first. However, nothing forbids us to still play with another representation of the 
particle tensor It turns out that this is the key for completing the construction 

which we summarize in the following theorem. 

Theorem 4. Every injective Matrix Product State is the ground state of a hopping 
theory which depends solely on the transfer matrix of the Matrix Product State. The 
Hamiltonian is the result of taking a long range and large system size limit of the parent 
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Hamiltonian 0 with correction matrix A set of m-particle eigenstates in the limit, 
are given by the symmetric super-positions, 

) •• • ) ) •• • ) ^rn}) 

PClSji Tl\... <C.Ti m 

summed over all permutations P. The basis states \(/)[Baj^,ni; na^] are given by, 

I 

Each tensor Ba^ is a function of Xa and ka, 

B(i) = - A^^Y^, F„ = (Id-e-'''r*)-i[X„] 

The energy of each m-particle eigenvectors is the sum of the energy of the individual 
particle, 

3 

The particle representations Xa and their respective energy are given in Lemma 

The states and respective energy of the theory are valid for a low-particle density 
which satisfies. 


0{D^\\\^) ^0,0 E (m — n + 1) 

\n=l 


L 

n 


N~^ /^2(n+l) 


^ 0, 0{m\N-^) 0 


Proof. Let us first fix m, and the particles representations and their 

respective momenta ki,...,km as given by Lemma (|^. One should first notice the 
following facts. First of all, one should see that one-particle representations with the 
tensor Ba^ is gauge-equivalent to Xa, 

\i>{Xa;ka}) = \tP{Bjf^;ka}) (A.ll) 

Secondly, the choice of gauge was chosen in such way that. 


Wki,k2,aY,AC 


jM 




Ebc 




'Al = o 




^3 


It is necessary for us to compute the overlap between these new basis states 

mBf.KY;.,), 


ni<...<nr7 


\(l){Bf = ... Ad-iSpAd+i ... Bff”^A^^+^ ... \i) 


With some algebra, or by using the gauge equivalence (A.ll), we derive for k = ka = kb, 

E (sf’Bf) = Tr 
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Using these simple properties of the tensor Ba ^ and particle representations Xa, one 
easily checks, 


m 


It follows that such basis states with different particle numbers are orthogonal, 






= o 


yYmin{mi,m2} 

mi+m2 

N 2 


Finally, this allows us the try and find some orthonormal basis in each subspace. We 
fix the particles Xa ^,..., Xa^ with respective momentum ki,, km and superimpose the 
basis states consisting of the permutations of the particles. 

; •• • ; 1^11 ••• ) ^m}) 

PCiSn 721 <...<77-771 

We verify that this choice yields the orthogonality, 

5 • • • 5 ^am 5 5 • • • 5 ^77i} |'0{^6i 5 • • • 5 ^hm 5 5 5 Q.m \) 

< ^ TT y I j{{qi-ki)np^^)^...+{qrn-krri)np^^^) 

— TYm il 

3 PeSm ^P(l)<--<?T'P(m) 

j ^ ^ 

The sought orthonormal basis has thus been found. 

One, now, needs to know the action of the Hamiltonian onto the m-particle states. 
Lemma’s 0 and 0 can be applied. 

The theorem is proven once we manage to show that the effective Hamiltonian of 
each m-particle subspace decouples. 


j 

Our choice of orthogonal basis indeed implies such result, 

5 ••• ? ^am"! ••• ? ^TTi} 5 ••• 5 ^bm"! 5 ••• 5 Q.m\) 

jl,j2 

With this the theorem has been shown. 


□ 
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We showed the symmetrized states could be considered as a type of eigenstates of 
the Hamiltonian. Unfortunately, this is only a small part of the spectrum. From our 
results in the previous lemma’s, we know that the n-particle spectrum is some complex 
superposition of the n-particle basis. However, in the case of normal unital quantum 
channels, this superposition turns out to be extremely simple. All n-particle energy 
levels are degenerate and a trivial basis can be found by simply taking the symmetric 
and anti-symmetric superposition of the n-particle basis. 

Corollary 5. When the transfer matrix is a normal unital quantum channel, in the 
long range and system size limit, the Hamiltonian 0 becomes exactly solvable with 
eigenstates 


) •• • ) ^a-m ) ) •• • ) ^m}) 

PClSji Til <C ... m 

with AiTa^yP) = iA(U). 


Appendix A.l. Normal unital trace-preserving Completely Positive Operators 


In general, we need not only information about the spectral properties of the transfer 
matrix, but also the algebraic relation between the eigenvectors. 

It turns out that normal unital trace preserving completely positive operators are 
an interesting class to study, as a lot of the computations, we had to endure for proving 
Theorem 0. simplifies. 


P = 


Tr(l) 


, r*[i] = 1, r[i] = 1 


It is important to stress that not any spectrum is allowed. Some, sometimes heavy, 
restrictions are imposed due to the complete positiveness. For example, if we assume 
that the eigenvectors have a particular group structure 

= Ug, Ug = (BjC— (A.12) 

This example is based on the lemma that spectral vectors around unit circle are of this 
form (|A.12|). 


27riKg 


Still, the following restrictions can be shown for T[Ug] = |Ag|e d JJg, 

27t [Kg - ga]' 


{!,...,21)}, 

9 

Writing, 


cos 


D 


> 0 


(A,13) 


r|.] = Ewn(c/t[,]) 


Equation (A. 13) follows for the necessary and sufficient for complete positiveness m. 

e[r] = El«>(/5l®r[|a)(/3ll>o 

a,(3 








Signatures of Lattice Excitations in Quantum Channels 


26 


This simple example illustrates that spectral properties of completely positives are not 
just arbitrary. There are physics restrictions. It is amazing that nonetheless, we can 
connect these with the spectrum of Hamiltonians, which are pretty much arbitrary. 

As presented at the end of the work, the appearance of bound states, is speculated 
from the breaking down of the convergence towards the one-particle subspace. 


lim ei{Xa,k) 

L—)-oo 


\H\tf{Xa,k})^2{L - ak,a)\HXa,k})\\l 


^k.a ff® 


1 - ed^-'^^a]|Ax, 


(A.14) 

(A.15) 


The following conclusions can be taken. If for all momenta k, ei(A’a, k) 0, then 
\ip{Xa,k}) are eigenstates with energy Re 1/(1 — 

In the case of bound states lying in the 2— and 1— particle, i.e. \'ip{Xa,k}) + 
Cgli^iXa, k — q; q}, we need a convergence of the 3-particle errors. 


e{Xa,q),e{Xi3,q) 0 
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